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Abstract 



A mass scale /x c representing the boundary between effective theories of strong interactions and perturbative 
QCD (infrared and ultraviolet regimes) exists when the /3-function is characterized by a simple pole. This 
behaviour, which is known to occur in N = 1 supersymmetric Yang-Mills theories, leads to an infrared attractor 
Oh! m the evolution of the coupling constant, with the mass scale fi c of the attractor providing a natural boundary 

between the infrared and ultraviolet regimes. It is demonstrated that [2|2], [3 1 1] and [1|3] Pade-approximant 
versions of the three- flavour (rif = 3) QCD /3-function each contain a simple pole corresponding to such an 
infrared attractor. All three approximants, separately considered, are seen to lead to nearly equivalent estimates 
for the mass scale fj, c . 



Although QCD is well-understood both qualitatively and quantitatively as a perturbative theory for the strong 
interactions, we have a surprisingly limited amount of information about the infrared boundary of its perturbative 
domain. Such a boundary (which is also anticipated by arguments in which hadrons and quarks are dual field 
variables for weak and strong QCD [Q) must exist to separate effective strong interaction theories (e.g., chiral 
perturbation theory, linear sigma model, etc.) from the higher-momentum region characterised by the perturbative 
quantum field theory of quarks and gluons. Consequently, one can argue that the infrared behaviour of the 
perturbative QCD coupling should not be characterised by smooth evolution to an infrared-stable fixed point, but 
rather by behaviour that would clearly separate the infrared domain of effective strong-interaction theories from 
perturbative QCD. Indeed, there exists both lattice j2|, analytical ||, and Pade-approximant j|, || corroboration 
for the idea that an infrared-stable fixed point does not occur within QCD unless the number of active fermion 
flavours contributing to the evolution of the QCD couplant is substantially larger than three. 

A clear demarcation between the infrared and ultraviolet regions would occur if the nj = 3 QCD /3-function were 
characterised not by a positive zero corresponding to an infrared-stable fixed point, but rather by a positive pole, 
an infrared-attractive singularity in the /3-function occurring at a momentum scale which necessarily corresponds 
to a lower bound on the domain of the running QCD coupling constant (assumed here to be real). Such behaviour, 
schematically presented in Fig. [A is known to characterise the exact /3-function for N — 1 supersymmetric Yang- 
Mills theory in the NSVZ renormalisation scheme ||, a "supersymmetric gluodynamics" whose /3-function pole 
forms an infrared-attractive terminal point for coupling-constant evolution within the theory's asymptotically-free 
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phase Jt| . There exists evidence from Pade-approximant methods that similar dynamics may also characterise the 
infrared region of QCD Q]. In the present note, we discuss whether such methods can provide information as to 
the actual mass scale which separates the infrared and ultraviolet domains of strong interaction physics, as well as 
the strong-interaction couplant magnitude characterising this mass scale. 

Weighted asymptotic Pade-approximant (WAPAP) methods have been utilised to estimate the n/-dependence 
of the five-loop contribution (^4) to the QCD /^-function, defined here for x(fi) = a s (/z)/7r to be 



H 2 j^ = -J2^x k+2 ^-(3ox 2 (l + R 1 x + R 2 x 2 + R 3 x 3 + R 4 x' l + ...) . (1) 



da; 

k=Q 

For rif — 3, the known MS /3- function coefficients in (|]) are 



9 3863 
£0 = 4 , 0i=4 , P2 = ^ , fo = 47.2280 . (2) 

The WAPAP estimate of the n /-dependence of the five-loop term (inclusive of quadratic-Casimir contributions to 
P3) is found to be ||] 

1 —3 
04 = jg (7.59 x 10 5 - 2.19 x 10 5 n/ + 2.05 x 10 4 n 2 f - 49.8^ - lMnj) 278 , (3) 

i?4 = ^^124 . (4) 

Po 

This estimate may be used to construct non-trivial N + M = 4 [iV|Af]-Pade approximants that reproduce the 
known coefficients Ri~R 3 and the estimate (^) for i? 4 within the [3- function series of ([!]): 

= 9 j2 [l- 5.4983a; - 1.9720a; 2 ] 
H y; 4 [1 - 7.2762a; + 6.4923a; 2 ] KJ 

m (x) = t 1 - 5 - 7397 "] (6) 

H V ; 4 [1- 7.5174a; + 8.8933a; 2 - 3.1896a; 3 ] V ; 

ram,, 9 Jl- 4.1155a; - 6.0058a; 2 - 5.3589a; 3 l 

- ~r ~ [1 - 5.8932,] ~ (7) 

All three approximants above exhibit a positive pole that precedes any positive zeroes, consistent with infrared 
dynamics analogous to those of NSVZ supersymmetric gluodynamics, as discussed above.^ This first positive 
pole, which corresponds to the couplant magnitude at the infrared-boundary momentum scale /z c , is surprisingly 
comparable for all three approximants: 

(8) 

The fact that all three approximants exhibit a positive pole of nearly equivalent magnitude (which precedes any 
positive zeroes) is indicative that such a pole is not likely to be an artefact defect pole jlTJ , but rather a manifestation 
of a true pole within the underlying all-orders /3-function. It is to be noted that the pole couplant magnitude, as 
estimated in (^|) , is itself sufficiently small for perturbative physics to remain viable near the infrared boundary. This 
behaviour is very different from an infrared-slavery scenario in which the QCD couplant grows non-perturbatively 
large in the vicinity of a deep-infrared Landau singularity 

Similar consistency of poles obtained from differing Pade approximants to the perturbative /3-function has 
already been shown to characterise supersymmetric gluodynamics in both NSVZ and in DRED renormalization 

1 The fact that a positive pole is always found to precede any positive zeroes for all three approximants has already been established 
M for arbitrary values of H4 . 
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= 0.160 
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X (fl c ) 
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[3|1] : 


X (fj, c ) 
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schemes [jO] . We reiterate that the pole in the former of these two schemes is known to occur from the all-orders 
/3-function expression, whether derived via NSVZ instanton calculus || or via imposition of the Adler-Bardeen 
theorem upon the supermultiplet structure of the theory [ p[ p"2| . 

One can utilise the infrared-attractive couplant values (JsJ) in order to obtain separate estimates of the infrared 
boundary /z c for each approximant considered. For specific n/ = 3 [JV|Af] approximant versions of the QCD 
/3-function, one finds that 



M fl M J=m T exp 



If dx 



2 J fi N \ M \(x) 



(9) 



We utilise the approximants within the integrand of (g), as well as the respective values (||) for [2|2], [1|3], 

and [3|1] approximant values of x(/i c ) for the upper bound of integration. For the lower bound of integration, 
we assume af s (m T ) = ttx (m T ) = 0.33 ± 0.02, consistent with recent analyses [CsJ . We then obtain via (g) the 
following values for the infrared-boundary momentum scale: 

^ 2|21 = 1.14 ± 0.11 GeV 

/4 1|3] = 1.13 ±0.11 GeV (10) 
/4 3|1] = 1.09 ±0.11 GeV 

These results not only exhibit remarkable consistency with each other, but also support the identification of 
perturbative QCD's infrared boundary with a momentum scale at (or somewhat above) the mass scale characterising 
nucleons and the vector meson octet. This picture is quite different from the usual one of an O(300 MeV) value for 
the Landau singularity Aqcd characterising coupling constant evolution via the truncated /3-function series. We 
reiterate that the estimates (||) and (|l^) for the infrared terminus of the couplant evolution within the asymptotically 
free phase of rif —3 QCD are consistent with both the applicability of controllably-perturbative QCD down to 
0(1 GeV) momentum scales, as well as the necessity for alternative descriptions (e.g. effective Lagrangians and 
hadronic field variables) to characterise sub-GeV (or sub-47r/„. JljJ) strong-interaction physics. 

The results (g) and @ are, of course, sensitive to input information. An across-the-board 150 MeV decrease 
from ( |io| ) in the estimated range of fi c is seen to occur if we choose the three-flavour threshold at /it = m c (m c ) = 
1.25 GeV @. More significantly, the set of values © for \i c relies ultimately on ref. J§'s WAPAP estimate (§,|) 
for the five-loop contribution to the QCD /3-function. The rif = 3 estimate obtained via (^) involves relatively 
small differences of large numbers. Consequently, the individual coefficients of powers of nj are likely to be more 
accurate than the estimate of @4 and i?4 obtained by evaluating the full polynomial in 

To test the uniformity of the infrared-boundary mass scale obtained via different approximants, let us choose 
i?4 for each approximant so as to ensure the occurrence of a positive pole at a s (/i c ) = 7r/4 (i.e.. at x = 1/4). This 
particular choice is motivated both as the critical value of the strong coupling for chiral symmetry breaking []l7| as 
well as by the couplant value characterising Nambu-Jona-Lasinio and linear-sigma-model approaches toward the 
generation of a dynamical quark mass |lq] . 

The appropriate set of approximants for the rif — 3 /3-function series 



2 d.X 9 n 



d^ 2 



1 + 16 x + 3863^ 2 + 2QQ902x 3 + RiX 4 
9 864 



(11) 



2 Such is the case for a prior asymptotic Pade-approximant prediction of the four-loop term [ [lq ] , /3^ red = {23600 — 6400n f + 350n 2 + 
1.499n^}/256, whose polynomial coefficients are in good term-by-term agreement with those of the subsequent exact calculation [^j, 
ptrue = {29243.0 - 6946. 30^ + 405.089n 2 , + 1.499n?}/256. When n f = 3, however, the estimated value /3f red = 29.65 differs from 
ptrue _ 47 228 by a relative-error (37%) whose magnitude is a factor of two or more larger than that of the relative error characterising 
each estimated polynomial coefficient (19%, 7.9%, and 14%, respectively). 
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a s (ra T ) = 0.31 


a s (m T ) = 0.33 


a s (m T ) = 0.35 




785 McV 


874 McV 


960 McV 


mL 1 ' 31 


778 McV 


867 McV 


952 McV 


mL 3111 


775 McV 


863 McV 


948 McV 



Table 1: Values of /x[ ; Ar ' M ', as obtained via (|9|), for [7V|M] approximants to the nj = 3 QCD /3-function. The 
(unknown) five-loop contribution i? 4 to the perturbative series ( |ll"| ) is chosen (see (|l5|)) to ensure a pole at x (fi c ) = 
1/4 for all three Pade-approximant cases. 



with i?4 arbitrary is 

9 2 [l + (7.19456 - 0.102610i? 4 ) x + (-11.3292 + 0.0756438i? 4 ) x 2 ] 



(3™(x) 



(3 [m {x) = - A 



4 [1 + (5.41678 - 0.102610i? 4 ) x + (-25.4301 + 0.258062i? 4 ) x 2 ] 

[1 + (5.80845 - 0.0933552i? 4 ) x] 



(12) 
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^\x) 



4 [1 + (4.03067 - 0.0933552i? 4 ) x + (-11.6367 + 0.165965i? 4 ) x 2 + (-18.3242 + 0.122349i? 4 ) x 3 } ' 

(13) 

9 2 [1 + (1.77778 - 0.0476412i? 4 ) x + (4.447106 - 0.0846954i? 4 ) x 2 + (20.9902 - 0.213007i? 4 ) a; 3 ] 



[1 - 0.0476412i? 4 a;] 



(14) 



The first positive pole for all three approximants is seen to precede any positive zeroes and is clearly dependent 
on the value of R4. One can easily show that all three approximants acquire a positive pole at x = 1/4 for very 
similar values of i? 4 : 



[2|2] 
[1|3] 
[3|1] 



R 4 = 80.307 
Ri = 89.925 
R A = 83.961 



(15) 



Given a s (m T ) and x (fi c ) = 1/4, as ensured by these respective choices for i? 4 , one can utilise (|) to estimate 
corresponding mass scales /i c for the three approximants. In Table [|, such values are obtained via couplant 
evolution from three different choices for a s (m T ) in the range a s (m T ) = 0.33 ± 0.02 JljJ. In all these estimates, 
three-flavour coupling constant evolution is assumed to be valid below m T . 

Table ^ shows striking uniformity in the values for the infrared boundary mass scale fi c obtained via three 
distinct Pade-approximants. Moreover, the three approximants generate values for \i c very near the p-meson mass 
when the lowest value for a s (m T ) is chosen. 

Indeed, Figure || demonstrates that the infrared boundary for all three approximants is nearly equivalent even 
if R4 is allowed to vary arbitrarily. The corresponding mass scale /j, c , as plotted in Figure is obtained for each 
approximant by substituting ([j^-pr] ) into the integrand of (^), with [|)'s upper bound of integration identified with 
the first positive poles of (p^-|l4|) .^For a given choice of i? 4 , Figure [2 shows that the infrared-boundary mass scales 
characterising all three approximants are surprisingly close. The O(600 MeV) lower bound evident from the figure 
for all three approximants is particularly striking. Different Pade approximants to the rif — 3 QCD /3-function thus 
appear to be quite consistent in the infrared dynamics they predict, suggestive that similar pole-driven dynamics 
may characterise the infrared boundary of QCD itself. 

The above results have all been obtained in the MS scheme, an explicitly gauge-invariant renormalization 
procedure for which the /3-function has been explicitly calculated to four-loop order, thereby facilitating our use of 

3 When R4 becomes negative, the [3 1 1] approximant version of the nj = 3 /3-function no longer exhibits a positive pole (or any 
positive zeroes, as would be the case with an infrared stable fixed point), but still exhibits a Landau singularity extractable from (B() 
provided the upper bound of integration in (H) is replaced by infinity. The [3|1] approximant curve in Figure displays the mass scale 
for this Landau singularity when R4 is negative. The figure shows this mass scale to be consistent with the mass scales characterising 
the simple poles within the [2|2]- and [l|3]-approximant versions of the /3-function for the same negative values of R4. 
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Pade approximants to the series ( [j]). Moreover, perturbative QCD phenomenology is more readily available (and 
more likely to be corroborated) in MS than in other renormalization schemes. The question remains, however, as to 
whether the infrared boundary we observe is a peculiarity of the MS scheme we employ. In eq. (Q), the coefficients 
i?2, i?3, -R4, . • • , (corresponding to (3k with k > 2) are all scheme-dependent, and therefore negotiable within the 
context of formal perturbative QCD. Indeed, a QCD renormalization scheme proposed by 't Hooft |L9| in which 
[3k — for k > 2 is guaranteed by construction to be free of /3-function polesj^J 

Recent work pOj has developed a procedure by which the leading renormalization-scheme dependence (i.e., 
explicit dependence on R2) can be eliminated from a perturbative QCD result. This approach leads ultimately to 
/3-functions characterized by the remaining renormalization-scheme parameters (i.e., R3, . . .). A reasonable set 
of Pade approximants to /3-functions extracted by this procedure is shown in to exhibit poles of a magnitude 
(0.29 < x < 0.33) sufficiently large for chiral symmetry breaking at the infrared boundary of QCD. 

Acknowledgments: We are grateful for useful discussions with V.A. Miransky, and for support from the Natural 
Sciences and Engineering Research Council of Canada (NSERC). 
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Figure 1 : Qualitative behaviour of the running coupling x (/i) in an asymptotically- free theory with an infrared 
attractor devolving from a simple pole in the /3-function. 
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Figure 2: Dependence of the infrared-boundary scale ^ c on the five- loop (3- function coefficient R4 — fli/fio, based 
upon n/ = 3 evolution of the couplant from an initial value of a s (m T ) = 0.33. We only plot values of R4 less than 
the WAPAP estimate (Q); if R4 is allowed to increase much past this value, the first positive poles of (|l2]-|l4|) are 
soon seen to be reached at values of /i exceeding the four- flavour threshold fit- 
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